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ABSTRACT:  Diffusion  of  adsorbates  on  transition  metal  nanoparticles  is  a  precursor  process  for 
heterogeneously  catalyzed  reactions,  and  as  a  result,  an  atomistic  understanding  of  the  diffusion 
mechanism  is  very  important.  We  systematically  studied  adsorption  and  diffusion  of  atomic  and 
diatomic  species  (H,  C,  N,  O,  CO,  and  NO)  on  nanometer-sized  Pt  and  Cu  nanoparticles  with  different 
sizes  and  shapes  using  density  functional  theory  calculations.  We  show  that  nanoparticles  bind  adsorbates 
more  strongly  than  the  corresponding  extended  single  crystal  metal  surfaces.  We  find  that  there  is  a 
Bronsted-Evans-Polanyi  type  linear  correlation  between  the  transition-state  energy  and  the  initial-state 
energy  for  adsorbate  diffusing  across  the  edges  of  Pt  and  Cu  nanoparticles.  We  further  show  that  the 
barrier  for  adsorbate  diffusion  across  the  nanoparticles  edges  can  be  estimated  by  the  binding  energy  of 
the  adsorbate  on  the  nanoparticles.  These  results  provide  useful  insights  for  understanding  diffusion- 
mediated  chemical  reactions  catalyzed  by  transition  metal  nanoparticles  which  are  widely  used  in 
heterogeneous  catalysis. 
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Diffusion  of  adsorbed  species  on  transition  metal  surfaces  is  an  important  process  for  thin-film  and 
nanostructure  growth  and  for  heterogeneous  catalysis,  among  others.'""^  In  connection  to  the  latter,  such 
surface  diffusion  is  the  precursor  process  for  the  subsequent  heterogeneously  catalyzed  reactions,  thereby 
affecting  catalytic  chemistry. 

Nilekar  et.al.  have  shown  that  the  transition  state  energy  for  the  diffusion  of  a  variety  of  catalytically 
relevant  adsorbed  species  on  single  crystal  transition  metal  surfaces  correlates  well  with  the  binding 
energy  of  the  adsorbate  on  the  respective  surface.^  This  correlation  is  a  Bronsted-Evans-Polanyi  (BEP)- 
type  of  correlation,  similar  to  other  BEP  correlations  established  earlier  for  surface-catalyzed  bond- 
breaking/bond-making  reactions. The  universal  BEP-type  correlation  is  independent  of  the  nature  of  the 
adsorbed  species  and  that  of  the  metal  surface.  Eor  a  certain  class  of  surface-catalyzed  reactions,  the 
existence  of  a  BEP-type  correlation  reflects  a  similarity  between  the  geometry  of  the  transition  state  and 
that  of  the  final  state,  resulting  in  a  linear  relationship  with  a  slope  near  to  unity.^  Eurthermore,  for  surface 
diffusion  on  transition  metal  surfaces,  a  rule  of  thumb  on  the  diffusion  barrier  was  found,  i.e.,  the 
diffusion  barrier  can  be  estimated  by  ca.  12%  of  the  magnitude  of  the  binding  energy  of  the  adsorbate.^ 

Compared  with  bulk  metals  and  extended  metal  surfaces,  nanometer-sized  metal  nanoparticles  have 
unique  physical  and  chemical  properties.  Eor  example,  though  bulk  Au  and  extended  Au  surfaces  are 
quite  inert,'®  very  small  Au  nanoparticles  are  active  for  CO  oxidation,""'^  water-gas-shift  reaction,'® 
formic  acid  decomposition,'^’'^  and  selective  oxidation  of  styrene  towards  its  epoxide. The  unique 
properties  of  small  metal  nanoparticles  have  been  generally  attributed  to  the  presence  of  under¬ 
coordinated  sites  (e.g.,  the  bridge  and  top  sites  on  the  edges,  the  vertex  site  at  the  corners),  possible 
quantum-sized  effects,  specific  shape  of  the  nanoparticles,  strain  developed  on  their  facets,  and  the 
interaction  between  the  nanoparticles  and  their  support.^®'^®  Eor  transition  metal  nanoparticles  catalyzed 
reactions,  adsorption  and  diffusion  of  species  on  the  nanoparticles  are  important  processes;  yet,  several 
related  questions  remain  open:  (1)  How  strong  do  the  nanometer-sized  transition  metal  nanoparticles  bind 
the  adsorbed  species  compared  to  the  corresponding  extended  single  crystal  metal  surfaces?  (2)  Eor 
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diffusion  across  the  edge  of  nanopartieles,  does  the  BEP-type  eorrelation  between,  for  instanee,  the 
transition-state  energy  and  the  initial-state  energy  hold?  (3)  Is  there  a  way  to  estimate  the  energy  barrier 
for  erossing  the  edge  of  a  nanoparticle,  where  the  edge  eonneets  two  adjaeent  faeets  of  the  nanopartiele, 
from  the  binding  energy  of  the  adsorbate  on  a  facet  of  the  nanopartiele?  The  answers  to  these  questions 
are  of  paramount  importanee  for  understanding  the  physics  and  chemistry  oeeurring  on  the  nanopartieles. 
However,  due  to  the  exeessive  computational  cost  of  nanopartiele  modelling  through  eleetronie  strueture 
theoretieal  methods,  a  quantitative  study  with  aeeurate  density  funetional  theory  (DFT)  ealeulations  is 
still  missing. 

Here,  we  perform  a  systematie  DFT  study  of  the  adsorption  and  diffusion  energetics  of  selected 
adsorbed  speeies  on  transition  metal  nanopartieles.  We  ehoose  four  atomie  speeies  and  two  diatomie 
moleeules  (H,  C,  N,  O,  CO,  and  NO)  as  adsorbates  and  study  the  adsorption  and  diffusion  of  these 
adsorbates  across  the  edges  of  Pt  and  Cu  nanopartieles  with  different  sizes  and  shapes:  1 .07  nm  Ptss  and 
1.62  nm  Pti47  euboetahedral  nanopartieles,  and  0.96  nm  Cuss  and  1.46  nm  Cui47  icosahedral 
nanopartieles.  We  show  that  nanopartieles  bind  adsorbates  more  strongly  than  the  corresponding  extended 
single  erystal  surfaees;  a  BEP-type  eorrelation  between  the  transition-state  energy  and  the  initial-state 
energy  holds  for  these  six  adsorbates  diffusing  aeross  the  edge  of  Pt  and  Cu  nanopartieles;  and  the 
diffusion  barrier  aeross  the  edge  of  Pt  and  Cu  nanopartieles  ean  be  estimated  by  the  binding  energy  of  the 
adsorbate  on  Pt  and  Cu  nanopartieles. 

All  ealeulations  were  performed  using  the  Vienna  ab-initio  Simulation  Paekage  (VASP)  eode^^’^^  based 
on  density  funetional  theory.  The  projeetor  augmented  wave  (PAW)  potentials^^’^®  were  used  for  eleetron- 
ion  interaetions  and  the  generalized  gradient  approximation  (GGA-PW91)^^  was  used  to  deseribe  the 
exehange-eorrelation  funetional.  The  eleetron  wave  function  was  expanded  using  plane  waves  with  an 
energy  cutoff  of  400  eV.  The  ealeulated  Pt  and  Cu  lattiee  eonstants  (3.99  and  3.64  A,  respeetively)  were 
used;  these  values  are  in  good  agreement  with  the  respeetive  experimental  values  (3.92  and  3.62  A, 
respeetively).^^  The  extended  Pt(lll)  and  Cu(lll)  surfaees  were  modeled  by  a  four-layer  slab  with  a 
(2x2)  unit  eell;  whereas  the  open  Pt(OOl)  and  Cu(OOl)  surfaees  were  modelled  by  a  six-layer  slab  in  a 
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(2x2)  unit  cell.  A  vacuum  layer  of  '-'12  A  was  used  to  separate  periodic  images  along  the  z  direction.  For 
Pt,  cuboctahedral  Ptss  and  Pti47  nanoparticles  were  used;  for  Cu,  considering  the  instability  of 
cuboctahedral  Cu  nanoparticles  up  to  2000  atoms,  icosahedral  Cuss  and  Cui47  nanoparticles  were 
selected.^^  The  Pt  and  Cu  nanoparticles  were  modeled  by  a  supercell  with  3-dimenisonal  periodic 
boundary  conditions  in  a  simple  cubic  cell.  A  large  vacuum  region  along  each  of  the  three  dimensions 
was  used  to  ensure  that  the  interactions  between  each  nanoparticle  and  its  periodic  image  are  negligible 
(the  distance  between  two  nearest  surface  atoms  in  neighbouring  images  is  at  least  16  A).  The  Brillouin 
zones  of  the  Pt(l  1 1),  Cu(l  1 1),  Pt(OOl),  and  Cu(OOl)  surfaces  were  sampled  using  a  (8x8x1)  k-point  mesh 
based  on  the  Monkhorst-Pack  scheme. The  Brillouin  zones  of  Pt  and  Cu  nanoparticles  were  sampled  by 
the  Gamma  point  only.  The  bottom-two  layers  of  Pt(lll),  Cu(lll),  Pt(OOl),  Cu(OOl)  were  fixed  during 
relaxation.  Ptss,  Pti47,  Cuss,  and  Cui47  nanoparticles  were  fully  relaxed  without  any  constraint.  All 
structures  were  fully  relaxed  until  the  Hellmann-Feynman  forces  acting  on  the  atoms  were  smaller  than 
0.02  eV/A.  The  climbing  image  nudged  elastic  band  (CI-NEB)  method^^  was  used  to  calculate  the 
diffusion  barriers  of  adsorbates.  The  binding  energy  (BE)  is  defined  as  BE=  Eads  -  Edean  -  Egas,  where  Eads, 
Eciean,  and  Egas  are  the  calculated  total  energies  of  the  slab  (nanoparticle)  with  adsorbate,  the  clean  slab 
(nanoparticle),  and  the  adsorbate  species  in  the  gas  phase,  respectively. 

We  first  calculated  the  binding  properties  of  H,  C,  N,  O,  CO,  and  NO  on  Ptss,  Pti47,  Cu55  and  Cui47 
nanoparticles.  On  cuboctahedral  Pt  nanoparticles  which  expose  both  (001)  and  (111)  facets,  we 
considered  different  high  symmetry  adsorption  sites  on  both  (001)  and  (111)  facets:  (a)  fee,  hep,  bridge 
(Bill)  sites  on  the  (111)  facet,  (b)  bridge  (BlOO),  top  (AlOO),  hollow  (hoi)  sites  on  the  (001)  facet,  (c) 
the  under-coordinated  bridge  (BE  10)  and  top  (AEIO)  sites  on  the  edge,  and  the  vertex  site  (AV)  at  the 
corner.  On  icosahedral  Cu  nanoparticles  which  expose  (111)  facets  only,  different  high  symmetry  sites  on 
the  (111)  facet  as  well  as  bridge  site  (BEl  1)  and  top  site  (AEl  1)  on  the  edge  between  two  adjacent  (111) 
facets,  and  the  vertex  site  at  the  corner  were  considered.  These  different  adsorption  sites  are  schematically 
shown  in  Eigure  1. 
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Figure  1.  (a)  Different  adsorption  sites  on  the  cuboctahedral  Ptw?  nanoparticle  (only  part  of  the  outermost  shell  of 
atoms  is  shown):  (1)  AlOO:  top  site  on  the  (001)  facet.  (2)  BlOO:  bridge  site  on  the  (001)  facet.  (3)  hol-side:  hollow 
site  (middle)  on  the  (001)  facet.  (4)  hoi:  hollow  site  (comer)  on  the  (001)  facet.  (5)  AEIO:  top  site  at  the  edge 
between  (001)  and  (111)  facets.  (6)  BElO-b:  bridge  site  (middle)  at  the  edge  between  (001)  and  (111)  facets.  (7) 
BElO-a:  bridge  site  (comer)  at  the  edge  between  (001)  and  (111)  facets.  (8)  AV:  vertex  site.  (9)  hcp-a:  hep  site 
(comer)  on  the  (1 1 1)  facet,  (a)  B1 1 1-a:  bridge  site  (close  to  corner)  on  the  (111)  facet,  (b)  fee  site  on  (1 1 1)  facet, 
(c)  Blll-b:  bridge  site  (middle)  on  the  (111)  facet,  (d)  hcp-b:  hep  site  (middle)  on  the  (111)  facet.  Red  and  blue 
spheres  denote  adsorbates  and  Pt  atoms.  Pt  atoms  on  one  (001)  facet,  one  (11 1)  facet,  and  the  edge  between  them 
are  highlighted  by  grey,  cyan,  and  pink  spheres,  respectively,  (b)  For  icosahederal  Cuhv  nanoparticle,  only 
adsorption  sites  on  the  (111)  facet,  on  the  edge,  and  at  the  vertex  are  available.  The  notations  of  BE  10  and  AEIO 
are  changed  to  BE  1 1  and  AE 1 1 ,  respectively. 


Summarized  in  Table  1  are  the  most  favorable  adsorption  sites  and  the  calculated  binding  energies  for 
H,  C,  N,  O,  CO,  and  NO  on  the  Pt  and  Cu  nanoparticles  examined  here.  For  comparison,  we  also 
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calculated  the  binding  energy  for  these  adsorbates  on  the  flat  Pt(lll),  Pt(OOl),  Cu(lll),  and  Cu(OOl) 
surfaces  and  the  respective  data  are  also  included  in  Table  1.  As  it  is  clearly  seen  from  Table  1,  on  the 
single  crystal  close-packed  (111)  surfaees,  the  fee  site  is  the  most  favorable  for  all  six  adsorbates  studied; 
whereas  on  the  open  single  crystal  (001)  surfaces,  the  bridge  and  hollow  sites  are  favorable.  Note  that  as  a 
general  trend,  the  open  (001)  surfaces  bind  adsorbates  more  strongly  than  the  elose-paeked  (111) 
surfaces,  primarily  beeause  metal  atoms  have  lower  coordination  numbers  on  the  (001)  surfaees. 

On  nanoparticles,  mueh  more  diversified  trends  were  found.  On  the  cuboctahedral  Pt  nanoparticles, 
and  for  H,  O,  CO,  and  NO,  the  under-coordinated  bridge  site  at  the  edge  (BE  10)  and  the  bridge  site  on  the 
(001)  facet  (BlOO)  are  very  stable.  For  H  on  Ptss  (denoted  as  H/Ptss — similar  notations  will  be  used 
thereafter),  the  most  stable  site  is  BlOO,  with  a  binding  energy  (BE)  of  -3.05  eV.  The  BE  10  site  is  slightly 
less  stable,  with  a  binding  energy  of  -2.94  eV,  identical  to  the  binding  energy  of  H  on  the  flat  Pt(OOl) 
surface.  For  O/Ptss,  the  BE  10  site  is  the  most  stable,  with  a  binding  energy  of  -4.75  eV,  followed  by  the 
hep  site  on  the  (111)  facet  with  a  binding  energy  of  -4.60  eV.  For  CO/Ptss,  the  binding  energy  on  BlOO, 
AV,  BEIO,  AEIO  is  -2.44,  -2.42,  -2.23,  -2.24  eV,  respeetively,  all  stronger  than  that  of  CO  on  the  flat 
Pt(OOl)  surface  (BE  =  -2.11  eV).  For  NO/Ptss,  the  binding  energy  on  BlOO,  BEIO,  and  AV  is  -2.61,  - 
2.56,  -2.56  eV,  respectively.  We  note  in  passing  that  the  BElO-b  site  (the  middle  bridge  site  on  the  edge) 
is  the  most  favorable  for  H/Pti47,  0/Pti47,  CO/Pti47,  and  NO/Pti47,  with  a  binding  energy  of  -3.04,  -4.66,  - 
2.32,  and  -2.54  eV,  respectively.  The  hep  site  near  the  eorner  is  the  most  stable  site  for  N/Ptss  and  N/Pti47, 
with  a  binding  energy  of  -5.37  and  -5.02  eV,  respectively.  For  C/Ptss  and  C/Pti47,  the  four-fold  hollow  site 
on  (001)  is  the  most  favorable,  with  a  binding  energy  of  -7.79  and  -7.90  eV,  respectively.  As  a 
comparison,  the  hollow  site  on  the  flat  Pt(OOl)  is  the  most  stable  adsorption  site  for  C,  with  a  binding 
energy  of  -7.71  eV. 
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Table  1,  The  most  favorable  binding  site  and  binding  energy  (BE,  in  eV)  for  H,  C,  N,  O,  CO,  NO  adsorbed  on  Ptss, 
Pti47,  Cu55,  Cui47  nanopaiticles.  Data  on  Pt(lll),  Pt(OOl),  Cu(lll),  and  Cu(OOl)  single  crystal  surfaces  are  also 
included  for  comparison. 


System 

Site!"] 

BE 

System 

Sitet^] 

BE 

H/Pt55 

BlOO 

-3.05 

H/Cu55 

fee 

-2.68 

C/Pt55 

hoi 

-7.79 

C/CU55 

4-foldt^] 

-5.66 

N/Pt55 

hep 

-5.37 

N/Cu55 

hep 

-4.32 

0/Pt55 

BEIO 

-4.75 

0/CU55 

hep 

-5.19 

CO/Pt55 

BlOO 

-2.44 

CO/CU55 

hep 

-1.19 

NO/Pt55 

BlOO 

-2.61 

NO/CU55 

hep 

-1.67 

H/Ptl47 

BElO-b 

-3.04 

H/Cu147 

hep-b 

-2.49 

C/Ptl47 

hoi 

-7.90 

C/CU147 

4-foldt'’] 

-5.26 

N/Ptl47 

hep-a 

-5.02 

N/Cu147 

hep-a 

-3.94 

0/Ptl47 

BElO-b 

-4.66 

0/Cui47 

hep-a 

-5.14 

CO/Ptl47 

BlOO-b 

-2.37 

CO/CU147 

AV 

-0.88 

NO/Ptl47 

BElO-b 

-2.54 

NO/CU147 

hep-a 

-1.27 

H/Pt(lll) 

fee 

-2.77 

H/Cu(lll) 

fee 

-2.53 

C/Pt(lll) 

fee 

-7.05 

C/Cu(lll) 

fee 

-4.83 

N/Pt(lll) 

fee 

-4.82 

N/Cu(lll) 

fee 

-3.73 

0/Pt(lll) 

fee 

-4.39 

0/Cu(lll) 

fee 

-4.79 

CO/Pt(lll) 

fee 

-1.76 

CO/Cu(lll) 

fee 

-0.92 

NO/Pt(lll) 

fee 

-1.96 

NO/Cu(lll) 

fee 

-1.22 

H/Pt(001) 

bri 

-2.94 

H/Cu(001) 

hoi 

-2.41 

C/Pt(001) 

hoi 

-7.71 

C/Cu(001) 

hoi 

-6.19 

N/Pt(001) 

hoi 

-4.58 

N/Cu(001) 

hoi 

-4.61 

O/Pt(001) 

bri 

-4.30 

O/Cu(001) 

hoi 

-5.11 

CO/Pt(001) 

bri 

-2.11 

C0/Cu(001) 

bri 

-0.82 

NO/Pt(001) 

bri 

-2.33 

N0/Cu(001) 

hoi 

-1.28 

[a]  See  Fig.  1  for  notation  of  different  adsorption  sites. 

[b]  Pseudo-four-fold  site  developed  from  the  distortion  of  the  local  (111)  facet. 
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On  the  icosahedral  Cu  nanoparticles,  we  found  that  the  hep  site  near  the  corner  (hcp-a  shown  in  Fig.  1) 
is  very  stable.  On  the  smaller  Cuss  nanoparticle,  the  most  favorable  site  is  the  hep  site  near  the  corner 
(hcp-a)  for  N,  O,  CO,  and  NO.  For  H/Cuss,  the  fee  site  is  0.13  eV  more  stable  than  the  hep  site. 
Interestingly,  for  C/Cuss,  we  found  that  C  placed  at  hep,  AEll,  and  BE  11  sites  as  initial  guess  to  the 
energy  minimization  procedure,  always  relaxes  to  pseudo-four-fold  hollow  sites  by  relaxing  the 
underneath  Cu  atoms  to  a  localized  square-like  Cu(OOl)  geometry — a  tendency  to  form  four-fold 
coordinated  C  as  on  the  flat  Cu(OOl)  surface.  On  the  larger  Cui47  nanoparticle,  H,  N,  O,  and  NO  all 
prefer  the  hep  sites.  Eor  CO/Cui47,  we  find  that  the  vertex  site  (AV)  and  edge  sites  (BEl  1  and  AEl  1)  are 
slightly  more  stable  than  the  hep  sites.  Similarly  to  C/Cuss,  for  C/Cui47,  the  pseudo-four-fold  hollow  site, 
which  developed  by  distorting  the  local  (111)  facet,  is  the  most  stable  site. 

By  comparing  the  binding  energy  of  different  adsorbates  on  Pt  and  Cu  nanoparticles  and  the 
corresponding  single  crystal  extended  surfaces  shown  in  Table  1,  we  found  that  nanoparticles  bind 
adsorbates  more  strongly.  Compared  with  the  single  crystal  close-packed  (1 1 1)  surface,  the  1.07  nm-sized 
cuboctahedral  Ptss  nanoparticle  binds  H,  C,  N,  O,  CO,  and  NO  stronger  by  0.28,  0.74,  0.55,  0.36,  0.68, 
and  0.65  eV,  respectively.  Compared  with  the  single  crystal  open  Pt(OOl)  surface,  the  Ptss  nanoparticle, 
binds  H,  C,  N,  O,  CO,  and  NO  more  strongly  by  0.11,  0.08,  0.79,  0.45,  0.33,  and  0.28  eV,  respectively. 
The  increase  in  binding  strength  on  nanoparticles  can  be  attributed  to  the  presence  of  under-coordinated 
sites  on  the  nanoparticles.  Indeed,  the  single  crystal  Pt(OOl)  surface  (with  a  coordination  number  of  8  for 
its  surface  metal  atoms)  binds  adsorbates  more  strongly  in  general  than  the  close-packed  Pt(lll)  surface 
(with  a  coordination  number  of  9  for  its  surface  metal  atoms).  On  Pt  nanoparticles,  the  edge  and  vertex 
sites  of  nanoparticles  have  even  smaller  coordination  numbers  (7  and  5  respectively),  and  thus  tend  to 
bind  adsorbates  more  strongly. 

On  the  larger  1.46-nm  sized  Pti47  nanoparticle,  the  binding  strength  increase  with  respect  to  single 
crystal  surfaces  is  less  pronounced  than  that  found  on  the  smaller  Ptss  nanoparticle.  Compared  with  the 
single  crystal  close-packed  (111)  surface,  Pti47 binds  H,  C,  N,  O,  CO,  and  NO  more  strongly  by  0.27, 
0.85,  0.19,  0.27,  0.61,  and  0.58  eV,  respectively.  With  the  binding  energy  on  single  crystal  Pt(OOl)  as  the 
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reference,  on  Pti47  the  binding  strength  of  H,  C,  N,  O,  CO,  and  NO  increases  by  0.10,  0.20,  0.44,  0.36, 
0.26,  and  0.21  eV,  respectively.  The  stronger  binding  on  the  smaller  nanoparticles  can  be  understood  by 
the  finite  size  effect.^^’^"^'^^  As  the  diameter  of  nanoparticles  increases,  their  properties  can  be  viewed  as 
converging  to  the  respective  single  crystal  limits. 

Similar  trends  are  found  for  adsorption  on  the  icosahedral  Cu  nanoparticles  which  expose  only  (111) 
facets.  Compared  to  the  single  crystal  close-packed  Cu(lll),  the  smaller  0.96-nm  sized  Cu55  particle 
binds  H,  C,  N,  O,  CO,  and  NO  more  strongly  by  0.15,  0.83,  0.59,  0.40,  0.27  and  0.45  eV,  respectively. 
The  larger  1.64-nm  sized  Cui47  nanoparticle,  the  binding  strength  increases  by  0.43,  0.21,  0.35  eV  for  C, 
N,  and  O  respectively;  whereas  the  binding  energy  changes  for  H,  CO,  and  NO  is  very  small.  We  note  in 
passing  that  the  flat  Cu(OOl)  surface  binds  C  more  strongly  than  both  Cu  nanoparticles.  This  is  because 
icosahedral  Cu  nanoparticles  do  not  contain  four-fold  hollow  sites  which  optimize  the  four-fold 
coordination  for  C.  The  stronger  binding  on  the  smaller  nanoparticles  than  on  the  larger  ones  can  be 
explained  by  the  particle  size  effect.^^’^"^'^^ 
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Having  calculated  the  binding  trends  of  these  adsorbates  on  Pt  and  Cu  nanopartieles,  we  now  study 
adsorbate  diffusion  aeross  the  edge  of  the  nanopartieles.  On  Pt  cuboetahedral  nanopartieles,  all  possible 
diffusion  paths  aeross  the  edge  from  the  (001)  faeet  to  the  neighbouring  (111)  faeet  were  investigated; 
while  on  Cu  ieosahedral  nanopartieles,  all  possible  diffusion  paths  aeross  the  edge  between  two  adjaeent 
(111)  faeets  were  analyzed.  After  eomparing  different  diffusion  paths  for  each  adsorbate,  the  minimum 
energy  path  (MEP)  aeross  the  edge  of  the  nanopartieles  was  determined  on  both  Pt  and  Cu  nanopartieles. 
Representative  examples  of  the  MEPs  for  NO  diffusion  aeross  the  edge  of  cuboetahedral  Ptn?  and 
ieosahedral  Cui47  nanopartieles  are  shown  in  Eigure  2.  The  rate-determining  step  (RDS)  of  the  diffusion 
aeross  the  edge  was  then  identified  as  the  step  along  the  diffusion  reaction  coordinate  eharacterized  by  the 
largest  aetivation  energy  barrier,  and  the  initial-state  energy  (BEis),  transition-state  energy  (BEts),  and 
edge  barrier  (Eedge,  defined  as  BEts  -  BEis)  were  obtained.  Eor  example,  for  NO  diffusion  across  the  edge 
of  Pti47  nanopartiele,  the  RDS  is  from  the  bridge  site  at  the  edge  between  (001)  and  (11 1)  faeets  (BElO-b) 
to  the  hep  site  (hep-b)  on  the  (111)  faeet,  with  BEis,  BEts  and  Ees  of  -2.54,  -2.01,  and  0.53  eV, 
respeetively.  The  eomplete  dataset  for  24  adsorbate-nanopartiele  pairs  is  tabulated  in  Table  2.  A  elose 
look  at  Table  2  reveals  that  almost  all  the  RDS  for  the  adsorbates  diffusion  starts  from  or  ends  at  the 
bridge  sites  at  the  edge  eonneeting  adjaeent  faeets  exeept  for:  (i)  H  on  Ptss  and  Pti47  where  the  RDS  ends 
with  the  top  site  on  the  edge,  (ii)  C/Cuss  where  the  RDS  is  from  one  pseudo-four-fold  hollow  site  derived 
from  the  hep  site  to  another  pseudo-four- fold  site  at  the  edge  derived  from  the  BE  11  site,  and  (iii)  O  on 
Cu55  and  Cui47  where  the  RDS  is  from  one  hep  site  to  another  hep  site  in  the  adjaeent  (11 1)  faeet. 
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Figure  2.  Potential  energy  surface  for  NO  diffusion  across  the  edge  of  cuboctahedral  Ptny  nanoparticle  (blue  line) 
and  icosahedral  Cum?  nanoparticle  (red  line).  The  atomic  strucures  are  schematically  shown  in  the  inset:  Ptu?  on 
the  right,  Cuuv  on  the  left.  Blue  sphere  is  N,  and  red  sphere  is  O  atom.  The  energy  is  reffenced  to  the  total  energies 
of  clean  nanoparticle  and  gas  phase  NO. 
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Table  2,  Binding  Energy  of  the  initial -state  (BEis)  and  transition-state  (BEjs)  for  diffusion  across  the  edge  of 
nanoparticles  for  different  adsorbates.  The  TS  corresponds  to  the  rate  determining  step  (RDS).  Eedge  gives  the 
activation  energy  barrier  for  that  diffusion  event.  All  energies  are  in  eV. 


System 

RDS 

BEis 

BEts 

Eedge 

H/Pt55 

BlOO  ->  AEIO 

-3.05 

-2.81 

0.24 

C/Pt55 

hep  hoi 

-7.75 

-6.72 

1.03 

N/Pt55 

hep  BEIO 

-5.37 

-4.77 

0.60 

0/Pt55 

BEIO  ->  BlOO 

-4.75 

-4.14 

0.61 

CO/Pt55 

B100->  AEIO 

-2.44 

-2.03 

0.41 

NO/Pt55 

BEIO  hol-side 

-2.56 

-1.87 

0.68 

H/Ptl47 

BElO-b  ^  AEIO 

-3.04 

-2.80 

0.24 

C/Ptl47 

hol-side  BElO-b 

-7.72 

-6.68 

1.04 

N/Ptl47 

hep-a-^  BElO-a 

-5.02 

-4.39 

0.63 

0/Ptl47 

BElO-b  hep-b 

-4.66 

-4.24 

0.42 

CO/Ptl47 

BElO-b  hep-b 

-2.32 

-1.97 

0.35 

NO/Ptl47 

BElO-b  hep-b 

-2.54 

-2.01 

0.53 

H/CU55 

hep  BEl  1 

-2.55 

-2.39 

0.16 

C/CU55 

4-fold  ^  4-fold-bI"] 

-5.66 

-5.19 

0.47 

N/CU55 

hop  BEl  1 

-4.32 

-3.92 

0.40 

0/CU55 

hop  hop' 

-5.19 

-4.60 

0.59 

CO/CU55 

hop  BEl  1 

-1.10 

-0.96 

0.14 

N0/CU55 

hop  BEl  1 

-1.67 

-1.36 

0.31 

H/Cu147 

hop-b  BEl  1-b 

-2.49 

-2.44 

0.05 

C/CU147 

hop-b  BEl  1-b 

-5.24 

-4.89 

0.35 

N/Cu147 

hop-b  BEl  1-b 

-3.92 

-3.79 

0.13 

0/Cui47 

hop-a  hop-a' 

-5.14 

-4.61 

0.53 

CO/CU147 

BE-b  hop-b 

-0.82 

-0.76 

0.06 

NO/Cu147 

hop-a  ->  BEl  1-a 

-1.27 

-1.18 

0.10 

[a]  Pseudo-four-fold  hollow  site  at  the  edge  close  to  the  comer,  [b]  hep  site  in  the  adjacent  (111)  facet. 
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BE,3(eV) 


Figure  3.  (a)  Binding  energy  of  the  transition  state  (BEts)  versus  binding  energy  of  the  initial  state  (BEis)  for  the 
rate-determing-step  of  adsorbate  diffusion  across  the  edge  of  the  Pt  and  Cu  nanoparticles.  Data  shown  with  empty 
symbols  are  for  Pt  nanoparticles  (Ptss  black,  Ptw?  red).  Data  shown  with  fdled  symbols  are  for  Cu  nanoparticles 
(Cu55  blue  and  Cum?  magenta).  The  best  fit  line  gives  BEts  =  -0.03  -l-  0.88  BEis  with  R^=0.99.  (b)  The  barrier  over 
the  edge  (Eedge)  plotted  against  the  binding  energy  of  the  initial  state  (BEis)  for  the  RDS  of  adsorbate  across  the 
edge.  The  linear  regression  gives  Eedge  =  -0.03  -  0. 12  BEis  with  R^=0.66. 
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We  next  examine  the  relation  between  the  transition-state  energy  and  the  initial  state  energy  of  the 
RDS  of  all  24  data  points.  Plotted  in  Figure  3(a)  is  the  transition-state  energy  BEts  as  a  funetion  of  the 
initial-state  energy  BEis.  Remarkably,  there  is  a  good  BEP-type  correlation  between  BEts  and  BEis:  the 
best  fit  line  gives  BEts  (eV)  =  -0.004  +  0.89  BEis  with  i?^=0.99.  Eor  comparison,  on  Pt(l  1 1)  and  Cu(l  1 1) 
surfaces,  the  correlation  between  BEts  and  BEis  is  BEts  =  -0.10  +  0.87  BEis  with  i?^=0.99.  This  result 
suggests  that  a  great  similarity  exists  between  the  energetics  of  diffusion  events  on  small  metal 
nanoparticles  and  on  single  crystal  surfaces,  leading  to  a  nearly  universal  description  capable  of 
accounting  for  surface  diffusion  on  both  classes  of  surfaces. 

Einally,  we  turn  our  attention  to  address  the  question  on  estimating  the  barrier  for  adsorbates  to  cross 
the  edge  on  nanoparticles  from  the  binding  energy  of  the  adsorbates.  Erom  this  present  work,  on  the 
single  crystal  Pt(lll)  and  Cu(lll)  surfaces,  the  correlation  between  the  diffusion  barrier  Ea  and  the 
binding  energy  of  the  adsorbate  BEis  is  Ea  =  -0.1 1  -  0.13  BEis  with  i?^=0.60,  which  is  close  to  the  best  fit 
line  reported  previously.^  Eor  adsorbate  diffusing  across  the  edge  of  nanoparticles,  we  plotted  the  edge 
barrier  Eedge  as  a  function  of  the  initial-state  energy  BEis  in  Eigure  3(b)  and  fitted  the  data  using  linear 
regression.  Though  scatter  around  the  best  fit  line  exists,  the  edge  barrier  Eedge  can  be  correlated  with  the 
binding  energy  BEis  of  the  adsorbates  via  Eedge  =  -0.003  -  0.11  BEis  with  i?^=0.62,  i.e.,  the  edge  barrier 
Eedge  can  be  roughly  estimated  as  equal  to  1 1%  of  the  initial-state  energy  BEis.  This  indicates  that  the  rule 
of  thumb  on  surface  diffusion  barrier  on  single  crystal  transition  metal  surfaces  can  be  generalized  to 
adsorbate  diffusion  across  the  edge  connecting  adjacent  facets  on  transition  metal  nanoparticles. 

In  summary,  we  presented  a  systematic  density  functional  theory  study  on  adsorbate  adsorption  on  two 
small  nanoparticles  of  Pt  and  Cu  and  diffusion  of  those  species  across  the  edge  of  these  nanoparticles.  We 
showed  that  on  Pt  and  Cu  nanoparticles,  adsorbates  (H,  C,  N,  O,  CO,  and  NO)  bind  more  strongly  than  on 
the  respective  single  crystal  surfaces,  mainly  because  of  the  presence  of  under-coordinated  sites  on  the 
nanoparticles.  Size  of  the  nanoparticles  was  also  shown  to  play  a  role  in  determining  adsorbates’  binding 
energies  on  the  nanoparticles.  Remarkably,  we  found  that  for  adsorbate  diffusing  across  an  edge  of  the 
nanoparticles,  a  universal  Bronsted-Evans-Plonyi-type  correlation  between  the  transition-state  energy  and 
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the  initial-state  energy  holds.  Furthermore,  we  showed  that  one  ean  estimate  the  diffusion  barrier  for 
adsorbates  erossing  a  Cu  or  Pt  nanopartiele’s  edge  by  '-1 1%  of  the  binding  energy  of  the  adsorbates  at  the 
initial  state  of  the  diffusion  rate-determining-step. 
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